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Abstract. The effects of intense electromagnetic fields on the decay of quasistation- 
ary states are investigated theoretically. We focus on the parameter regime of strong 
laser fields and nonlinear effects where an essentially nonperturbative description is 
required. Our approach is based on the imaginary time method previously introduced 
in the theory of strong-field ionization. Spectra and total decay rates are presented for 
a test case and the results are compared with exact numerical calculations. The po- 
tential of this method is confirmed by good quantitative agreement with the numerical 
results. 
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1. Introduction 

Since tlie 1960s wlien lasers became a worldwide-used laboratory equipment, 
experimental and theoretical studies of nonlinear phenomena emerging in the interaction 
of intense electromagnetic fields with matter have seen continuous progress. One of 
the key phenomena in the physics of strong-field laser-matter interaction is nonlinear 
ionization. Although observed for the first time in 1965 [1], it continues to attract 
attention as a fundamental process that encodes pivotal information on quantum 
dynamics in the presence of intense electromagnetic fields. Another reason why nonlinear 
ionization remains in the scope of theoretical and experimental research is that it is 
a building block for numerous complex laser-induced and laser-assisted phenomena, 
including the generation of high order harmonics, nonsequential multiple ionization and 
a variety of effects in laser plasmas. 

Typically, the quantum system subject to an intense laser pulse is assumed to be 
stable in the absence of the field. In this sense, ionization is viewed as a laser-induced 
effect. This indeed holds for a variety of problems related to the interaction of laser 
pulses with atomic and molecular gases. Within this common picture, ionization differs 
qualitatively from laser-assisted phenomena, i.e., effects that also take place without the 
presence of an external field but can yet be essentially modified by it. Bremsstrahlung 
in external electromagnetic fields is such an example |2l [3]. Ionization can also be 
considered a laser-assisted effect under specific conditions, such as when the system is 
prepared in an excited state with a positive total energy (e.g., autoionizing states) or is 
under the action of a constant or slowly varying electric field and is thus free to decay 
by itself. In this case, an additional external field can influence the decay by modifying 
the electron spectra and the total probability. One may thus interpret such a process 
as being a laser-assisted decay (or, equivalently, ionization) from a quasistationary (QS) 
state. 

A QS state in the presence of an AC fleld can be in many cases formally treated 
as a true bound state affected by two flelds out of which one is static or quasistatic. 
From this point of view, laser-assisted decay is nothing but ionization in a two-color fleld. 
Nonlinear ionization and other related processes in the fleld of two or more superimposed 
intense laser pulses of different frequencies have been considered in theoretical studies. In 
particular, the interaction of laser radiation with negative ions in the presence of a static 
electric fleld was treated in detail in Ref . [1] , where the authors used the QS quasienergy 
states formalism, which combines the QS state and Floquet approaches. In addition, a 
number of theoretical works have considered two-color strong-fleld ionization. However, 
there is an essential difference in the statement of the problem that limits the possible 
application of the aforementioned results to the description of laser-assisted tunneling 
of a QS state: in the case of two-color laser-induced ionization, if one of the two flelds 
is turned off, the photoelectron spectrum remains an above-threshold-ionization (ATI) 
spectrum that consists of several ATI maxima separated by the photon energy. This 
holds for any ac fleld, irrespective of its frequency. If the second fleld is static (as 
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considered in Ref. then there is no spectrum along the direction of the electric 
field vector, since the respective momentum component is not conserved. In this case, 
the only observable of the problem is the rate of ionization, or the width of the QS 
state. When, however a true QS state is considered, in the absence of the laser field the 
spectrum has the form of a single narrow line centered at the real part of the QS state 
energy. None of the two aforementioned cases possesses this limit. 

The aim of this work is to develop a theoretical description of laser-assisted decay 
of QS states. For laser-induced ionization, a simple and simultaneously very efficient 
nonperturbative theory was proposed by Keldysh [5]. This approach, meanwhile known 
in its different realizations as the Keldysh- Faisal-Reiss model [HI [7] or the Strong-Field 
Approximation (SFA) [7] (for the present status of the SFA and its implementations 
see Refs. [H [9]), is a well-established tool in strong- field physics. The essence of 
the SFA is that it approximates the electron continuum by the Volkov states, which 
are exact solutions of the Schrodinger (Dirac) equation in the presence of a plane 
electromagnetic wave [IDl [H]. The effect of the binding potential on the electron 
continuum is simultaneously neglected. This key assumption of the SFA is intuitively 
natural and provides a way to a simple description of strong-field laser-atom interactions 
that in many cases is fully analytical. 

Here we adopt the general idea of SFA and modify it for the description of laser- 
assisted decay. For this task, we use the formulation of the SFA in terms of complex 
classical trajectories known as the Imaginary Time Method (ITM). This method 
was introduced in the early days of strong-field physics [13] to give a physically 
transparent formulation of the Keldysh method. In addition, the ITM also provides an 
efficient way to consider significant effects that the Keldysh model misses in its original 
formulatiorQ. The idea of using the ITM for the description of the tunneling of free 
particles through a potential barrier in the presence of an oscillating electric field was 
first proposed by Ivlev and Melnikov [19]. Here, using the ITM, we develop a more 
general and accessible theoretical description of the laser-assisted decay of QS in the 
presence of intense electromagnetic fields. Our approach recovers specific results of Ref. 
[T9] and goes beyond by providing not only qualitative but also quantitative tunneling 
probabilities that agree with exact numerical calculations. In the form presented here, 
restricted to the nonrelativistic limit, our method should be equally applicable for 
descriptions of different physical systems under standard semiclassical conditions. 

The general question about how AC electric fields can infiuence the decay of QS, i.e, 
tunneling, has been theoretically addressed in other fields. Most of the work we are aware 
of relates to the tunneling between superconductor films or micro- and nanocontacts 

X In particular, the ITM was used to account for the effect of fong-range Coulomb interaction on 
the ionization dynamics. It is currently proven in experiments and numerical simulations that the 
Coulomb interaction strongly affects nonlinear ionization. Coulomb corrections to the ionization 
amplitude derived within the ITM provide a quantitative description of total ionization rates \12\ 114] 
and photoelectron momentum distributions |15[ 1161 117] . The ITM and its applications were recently 
reviewed in Ref. [18]. 
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affected by a microwave, infrared or optical laser field. For such systems the effect 
is commonly referred to as photon-assisted transport (see |20] and references therein). 
Typically, in experiments on photon- assisted transport, low and moderate intensity laser 
fields are applied so that, although the coupling is nonlinear and multiphoton processes 
are present or even dominate, the characteristic number of absorbed or emitted photons 
is not very large. For such a moderately nonlinear regime of interaction, the Floquet 
method is known to be an efficient approach, and is routinely used [20]. However, when 
a large number of photons is involved, the Floquet approach becomes less efficient and 
more cumbersome numerically. On the other hand, the ITM is particularly relevant in 
the semiclassical domain where a typical number of photons is at least on the order of 
10 and can easily reach hundreds and even thousands. Thus, the approach we develop 
here can be viewed as a complementary one with respect to the Floquet method and 
may provide a new efficient computational tool in theoretical studies of photon-assisted 
transport in the highly nonlinear coupling regime. 

Another example of laser-assisted decay is beta-decay in the presence of 
electromagnetic fields |21]. The decay of elementary particles can be modeled within 
the picture of QS states only on a qualitative level. In the original work of Nikishov 
and Ritus [21], a more rigorous approach was used, based on the low-energy limit of 
the weak interaction theory modified to account for the effect of the electromagnetic 
field. To this end, the muon and the electron plane waves were replaced by Volkov 
waves. As will be shown later in this work, our approach reproduces an important 
qualitative conclusion of earlier studies [211 1221 1231 El], namely that depending on the 
parameters, two different regimes of decay are realized (i) when the spectrum is strongly 
affected without a modification of the total decay rate and (ii) when the rate of decay 
is also affected. These two regimes are also referred to as "exclusive" and "inclusive", 
respectively. Although Refs. [211 1221 (221 [21] consider the influence of laser flelds on the 
decay of elementary particles, the qualitative classiflcation of the interaction regimes 
holds for a variety of laser-assisted processes. 

The physical difference between the two regimes becomes clear if we notice that for a 
strong modiflcation of the spectrum no high-fleld intensity is actually needed, but only a 
large quiver (ponderomotive) energy. The latter can be achieved at low laser frequencies. 
A relatively weak but low-frequency laser fleld strongly affects the kinematics of the 
charged particle, accelerating or decelerating it after decay. This changes the flnal 
energy at the detector with almost no effect on the total decay probability. In contrast, 
in the "inclusive" regime, the particle's dynamics on the short-time scale corresponding 
to the subbarrier motion is also influenced, modifying the total probability. This indeed 
requires high laser fleld strengths. The particular expression of the critical fleld that 
delimits the "inclusive" regime depends on the investigated system. For beta-decay such 
flelds lie far beyond present experimental capabilities (see, e.g., Ref. [21] and references 
therein) and are probably in principle unachievable [25], whereas for atomic and solid 
state systems discussed in this work, both regimes are accessible. 

This paper is organized as follows. In Section |2] we introduce the SFA and adapt 
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it to describe the laser-assisted decay of QS states. The reformulation of SFA in terms 
of the ITM is used to derive the transition amplitude. With the help of the latter, we 
determine in section [3] the ATI-like spectra of a model one-dimensional (ID) problem 
where the QS state can tunnel through a rectangular barrier. The spectra are analyzed 
and compared with the results of the exact numerical solution of the time-dependent 
Schrodinger equation (TDSE). We conclude with a short summary and outlook. Atomic 
units h = TUe = \e\ = 1 are used throughout the paper. 



2. Basic equations 



2.1. Strong-field approximation (SFA) matrix element for quasistationary (QS) states 

We start from a short summary of the SFA that describes ionization from true bound 
states. Within the SFA, the transition amplitude between an atomic bound state |\l/o) 
of binding energy Eq = —I and a continuum state |\l/p) with an asymptotic momentum 
p is given by 



Msfa(p) = -^ J {^p\V{t)\<Ilo)dt, (1) 

— oo 

where the final state is approximated by the Volkov function, 

^p(r,t) = ^^^exp |zvp(t) -r-t J ep{t')dt'^ (2) 

and 

^pW = Vp(t)/2, with Vp(t) = p + A{t)/c 

are the electron time-dependent kinetic energy and velocity in the electromagnetic 
field, respectively, described by the vector potential A{t). V{t) is the interaction operator 
of the electron with the field of the electromagnetic wave and c the speed of light. In 
the dipole approximation, the electric field S{t) = —dfA/c and the vector potential 
depend only on time. It is convenient to simplify the notation by using the field-induced 
momentum rather than the vector potential, Pi?(t) = A(t)/c. 
Amplitude (II]) is relevant under the semiclassical conditions: 

i^o = - » 1 , = 1^ « 1 , (3) 

where Eq is the electric field amplitude, £ch = (2J)^/^ is the characteristic atomic field 
(for the ground state of hydrogen £ch = mle^ /h'^ = 5.14 x lO^V/cm) and cu is the laser 
frequency. The first strong inequality in ([3]) shows that the minimal number of photons 
required for ionization is large, hence the coupling is essentially nonlinear. The second 
inequality guarantees that the spatial scale on which the ionization amplitude forms 
is large in comparison with the atomic size (see a detailed discussion in Refs.[8], [T5]). 
Another frequently used dimensionless combination known as the Reiss parameter [7] is 
proportional to the ratio of the ponderomotive energy Up = (p|.(t))T/2 (where T is the 
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optical period) to the photon energy. For the hnearly polarized monochromatic field 
this reads 



Zf 



P C-o 



SF^Kl . (4) 



For ionization of atoms and positive ions by intense infrared and optical lasers, the 
conditions ([3]) are usually well satisfied and zp ^ 1. The integrand in ([T]) is then 
a rapidly oscillating function of time. This allows us to evaluate the integral by the 
saddle-point method, so that the amplitude can be written as a sum of contributions 
from all relevant stationary points to(p); 

exp (-«S'o(p,toa)) 



Mspa(p) = v^^E^(P^U 7/ (5) 

/(9/^o(P,^Oa) 



where S is the classical action 

+ 00 

^o(p,t)= / {^l{t')/2 + l}dt' (6) 



t 



and the pre-exponential factor V is the spatial matrix element of the interaction operator 
V . The saddle-point equation is of the form 

ai5'o(p,to) = v5(to)/2 + / = 0, (7) 

which shows that a saddle point is always complex for / > 0. The differential ionization 
rate is given by the squared modulus of ([5]). 

The fact that all roots of Eq. ([7]) are complex refiects the quantum nature of 
ionization. Consequently, phase 5*0 (p, to) in ([5]) is a complex quantity with a numerically 
large and negative imaginary part, hence under conditions ([3]) the ionization probability 
appears to be a highly nonlinear function of the laser field strength. The form of this 
nonlinear dependence is quantified by the value of the Keldysh parameter [5] 

y/2iu _ 1 



which is the ratio of the characteristic atomic momentum = yll to the field 
momentum -pp = Sq/u. Since the final state is approximated by the plane wave ([2]), 
the prefactor "P in ([5]) can be expressed via the Fourier transform of the bound state 
atomic wave function 00 (r)- For practical calculations, this means that the prefactor is 
a weak function of the final momentum and the field and atomic parameters compared 
with the highly nonlinear exponential, so that one may safely replace it by a constant 
as long as photoelectron momentum distributions and not the total ionization rates, are 
considered. Moreover, in the above formulation of the SFA, the so-called plain SFA, 
the expression of the prefactor is anyway incorrect, except for the case of ionization 
from a short-range well. The simplest form of the SFA transition amplitude is therefore 
given by Eq. (I5j) with V = const. On the qualitative level this rough approximation is 
in many cases sufficient. Strictly speaking, the SFA provides a quantitatively correct 
description of nonlinear ionization only for the exceptional case of a particle bound by 
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a zero-range potential. For short-range potentials, it is still a good approximation if 
the interaction operator V{t) is taken in the length gauge [26]. For atoms, where the 
electron-core interaction potential always presents a long-range Coulomb tail, the SFA 
prefactor V is essentially wrong in any gauge. In this case, to calculate it correctly and 
bring the SFA back to the quantitative level of description, the technique of Coulomb 
corrections was developed. For further details we refer the reader to Refs. [HI [151 127] . 

Let us now turn to the case when not a bound but a QS state is subject to an 
intense laser pulse. One would then expect the appearance of an ATI-like photoelectron 
spectrum with the significant difference that now the initial state energy Eq = > 

so that there is no gap between the initial state and the continuum. As a result, laser 
photons can also be emitted, not only absorbed, and the net number of absorbed photons 
can also be zero. Figure [T] sketches this qualitative difference between photoelectron 
spectra for stationary and QS states. Nevertheless, in strong fields we expect that the 
typical number of photons involved in the interaction is anyway large, hence SFA-like 
approaches should be suitable also for the description of ionization from QS states. With 
this assumption, we can introduce the amplitude of laser-assisted decay replacing the 
bound state wave function |\l/o) in (11]) by the QS state (Gamow's wave function) 

2 

vl/o(r, t) = 0o(r)e-*^*, E = E, + lE'^ = - iT/2 = ^- - tpop'o . 

Here Eq is the real part of the complex energy and F is the width that determines 
the decay rate. Following a common width limitation in the theory of QS states, we 
consider T Eq. If an even stronger limitation is satisfied and the width is small 
compared with all other characteristic frequencies of the problem, we may disregard the 
factor exp(— rt/2) in the integrand. Then the SFA ionization amplitude differs from the 
one for the true bound state in the spatial wave function of the initial state 0o(r) and 
by the fact that the initial state energy Eq lies in the continuum. 

This straightforward application of the SFA leads, however, to some difficulties, 
namely: 

(i) The spatial matrix element is divergent due to the exponential divergence of the 
spatial wave function of the QS state, 

(/.o(r) ~ e'fo^+Po'^ ^oo, r^oo, p'q ^ V/2pQ . 

This asymptotic exponential divergence at large distances well known in the theory 
of QS [281 129], originates from the approximate treatment of the decaying state as 
a stationary state and was noted in the pioneering work of Gamow |30] . 

(ii) Even if the phase is large and the saddle point method is applicable, the saddle 
points are real (since Eq > Q and we omit V <^ Eq ixi the saddle-point equation) 
and the stationary phase is also real. The field parameters then enter the tunneling 
probability only via the pre-exponential factor. In other words, the transition 
amplitude does not demonstrate the nonlinear dependence on the laser field strength 
and frequency typical for strong field phenomena. This means that even in a 
very weak field the probability of detecting the outgoing particle with an energy 




Figure 1. Sketch of ID short range potential U{x) including the potential well 
(0 < a; < a) and rectangular barrier (I). A smooth boundary to the rectangular barrier, 
including the triangular barrier of width A (II), is shown by the dashed-dotted line. 
Qualitative photoelectron spectra for strong-field ionization from a true bound state 
(B) and a QS state (QS) are presented. Dashed line: the infinite barrier used to obtain 
numerically with the help of the Schrodinger equation the ground-state wavefunction 
of the particle inside the well. For such a barrier with 6 oo. only the ATI spectrum 
(B) is present. See Section l3^ for further details. 



considerably different from Eq is not small, i.e. the field-free tunneling exponent 
does not emerge in the limit So ^ 0. Obviously, such a conclusion cannot be 
correct. 

Although the exponential divergence of QS state wavefunctions is itself not 
surprising and follows from the definition of a QS state, for the calculation of norms and 
matrix elements containing these divergent factors a regularization method is needed. 
Such a method was first proposed by Zeldovich [31] . However, for our purposes we will 
not use any regularization but apply instead another method for reconstruction of the 
correct prefactor, as explained in the following. 

The origin of the second difficulty becomes clear if we consider the structure of the 
continuum for the system shown in Figure [H For simplicity, in the following we refer to 
a ID system. At large distances the eigenfunctions are superpositions of the incoming 
and the outgoing plane waves: 

M^) = e"'^" + i^(p)e*^", x-^oo, E = 
while inside the well, x < a, 



(j).p{x) ~ A{p) sin 1^ J v{p, x')dx' j , v{p,x) = yp^ — 2U{x) . 
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For the momenta in the narrow vicinity of the QS state, p ^ the absolute value of 
the coefficient A{p) depends strongly on p having a sharp maximum at p = Pq. For 
eigenstates whose energy is sufficiently different from Eq, 

\E-Eo\:$>T, 0<E<Uo, (9) 

the wave function inside the well is exponentially small, 

b(p) 

A{p) e'^°^P\ Wo{p)= J v{p,x)dx. (10) 

a(p) 

Here a{p) and b{p) are the turning points of classical motion, hence the action Wq taken 
over the classically forbidden region is a purely imaginary value and iWo < 0. The 
coefficient A{p) is the semiclassical probability amplitude (calculated with exponential 
accuracy) of the particle to tunnel through the barrier formed by the potential U{x). 
Thus, for continuum states satisfying ([9]) the correct spatial matrix element should 
contain an exponentially small factor (ITOl) . whereas the SFA matrix element calculated 
with a plane wave final state function does not present this feature. We come to the 
conclusion that the problem with the SFA applied to QS states is that its plane wave final 
state Volkov function differs from the correct continuum wave function exponentially, 
exactly in that part of the position space that contributes most to the spatial matrix 
element. 

Taking this into account we may formulate how one should modify the SFA matrix 
element to make it appropriate also for the description of LAT from QS states; the 
spatial matrix element should be replaced according to 

{(j)p\V{x, t)\<f)o)/^JdiSo{p,toa) ^ A{vp)V{vp, t) , (11) 

where A{vp) is given by (ITOl) . Vp is the instant velocity Vp(t) = p + prif) and the new 
prefactor V{vp,t) will be defined below. Then, instead of ([5]) we obtain 

Msfa{p) = E ^(^'p. toa) exp [iWoiVp) - lSo{p, toa)] • (12) 

a 

The action Wq{vp) is complex and describes the field-free tunneling through the 
potential barrier U{x). The action So{p,to) is real (just as the saddle point to) and 
therefore accounts for the effect of the laser field on the particle after the tunneling. 
Thus, in this approximation the laser field only changes the particle's energy on its way 
from the tunneling exit to a detector. In other words, 5*0 accounts for the kinematic 
effect of the laser field that redistributes the particles in the energy space not affecting 
the total decay probability. This corresponds to the "exclusive" regime of interaction 
[2H [22| [23| [21]. At first sight, the tunneling exponent is affected by the laser field via 
the fact that the field-free action Wq is now taken at the instant velocity Vp at the 
saddle point. However, according to the saddle point equation ([71) Vp{to) = po, so that 
Wo{vp) = Wo{po) and the laser field dependence vanishes. 

We can now translate this formal description of the matrix element ( [T2l) into a 
simple physical picture which would allow us to determine the correct prefactor V{p, t). 
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dw{to) = Rodto = i?c 



dp, (15) 



The particle tunnels through the potential barrier the same way as it would without the 
laser field. At some time instant to, it emerges in the classically allowed domain having 
the initial velocity f (to) = Po and starts its motion in the laser fielcl§|. Then its final 
energy is given by 

E = py2 = {po-pF{to))V2, (13) 

so that each initial time to corresponds to a certain final energy. The inverse function 
to{E) is not necessarily single- valued. In the linearly polarized monochromatic field with 

£ (t) = So cos ut , p pit) = —pp smut , pf = £o/^ (14) 

we have -Emax = (po + VfY 1'^ and -Emin = (po — Pf)^/2 or zero; the spectrum consists 
of ATI-like peaks between the classical boundaries (CB) -Emm ^ E < E^^^^. The 
magnitudes of the ATI-like peaks vary slowly with the energy via the prefactor in 
Eq. f|T2|) . Outside the CB, the spectrum vanishes abruptly. 

Within the picture described by Eq. ( |T2l) . penetration of the particle through the 
well and its subsequent evolution are independent. Then the probability to tunnel out 
during the time interval dto is given by 

dto 
dp 

where Ro = Vq exp(— 2ImPVo) is the field-free tunneling rate and the derivative dp/ dto = 
—ppif) is calculated from (fT3|) . Note that we assume here that the tunneling probability 
can be written as a continuous function of the particle asymptotic momentum p. In 
reality, the spectrum consists of a discrete comb of peaks corresponding to absorption 
or emission of an integer number of laser photons. Our assumption thus implies that 
the characteristic number of peaks L in the spectrum is large, L ^ 1. According to 
Eq. (fT3|) . there are two limiting cases defined by the ratio between the field-free electron 
asymptotic momentum po and the field momentum pp. If Pf ^ Po then the classical 
boundaries of the spectrum are approximately Eo ±poPF so that the number of peaks 
is L ~ poPf/^- For stronger fields where Pf ^ Po the energy scale is determined by the 
ponderomotive energy and the number of peaks is of the order of the Reiss parameter 
(|4]). As will be discussed in the next section, for field values for which L is not much 
greater than unity, the approximation of the continuous spectrum contradicts energy 
conservation requirements. As a result, the total rate can only be calculated with some 
numerical error. 

The distribution ( IT5|) is divergent at the CBs p = Pmax/min where we have dp/ dto = 0. 
Such an integrable divergence near the CBs is typical for SFA-based descriptions, 
occurring at the final momenta for which the saddle-point method does not work due 
to cancellation of the second derivative of the action. This does not affect the total 
probability but renders the momentum distribution incorrect in the vicinity of CBs. To 
avoid this problem, the term that is proportional to the third derivative of the action has 
to be accounted for in the phase decomposition near the saddle point. The divergence 

§ Here we assume that the potential well is a short-range one. 
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is then replaced by a local maximum of the spectrum at the CB [22] with a relative 
height of the order of z]/^, where zp is the Reiss parameter In the simplest form, 
this regularization procedure reduces to the replacement 



dto 



f3 = £, 



1/3 



(16) 



dp ■ ^{dp/dtoy + /32 

Next, we take into account that usually there is more than one solution to the saddle- 
point equation, so that several toa? « = 1, 2, ... correspond to the same final energy. This 
leads to a coherent sum over all saddle-point solutions. The momentum distribution 
takes then the form: 

Vo exp [iWo{vp) - iSo{p, too)] 



dw{p) = \M{p)\'dp, M{p) = J2 



t=tOa 



(17) 



'dp/dto + i(3 

Clearly, this result misses two effects: (i) the influence of the laser field on the subbarrier 
motion is not accounted for and (ii) in the classically allowed domain, the effect of the 
potential U{x) is disregarded. The former effect becomes more and more significant 
when the laser intensity grows, whereas the latter is particularly important for potentials 
with a long-range tail, e.g. for the Coulomb potential. In the next subsection, we 
reformulate the amplitude ( IT7|) using the ITM and show that this new formulation 
provides a straightforward way to account for the two missing effects. 



2.2. Imaginary time method for quasistationary (QS) states 

The amplitude in Eq. fll7p can be equivalently reformulated in terms of classical complex 
trajectories. According to the ITM [18], a trajectory Xo(t) can be found along which the 
particle starts its motion at the complex time instant t = tg inside the well, x{ts) = 0, 
having the energy E = v'^{ts)/2 = Eq and arrives aX x = h when t = t^. Here h = 6(-E'o) 
is the outer classical turning point, U{b) = Eq. The trajectory satisfies the Newton 
equation 

x = v = —dU/dx. (18) 

The exit point x = b is separated from the well by the classically forbidden region; 
hence, the solution of Eq. (fT8|) satisfying the assigned initial conditions only exists in 
complex time, t = to + it. The action Wq in (fTOl) and (fTTl) can be represented as 

to 

WM = l{C + Eo)dt-pob, (19) 

ts 

where C = — U{x) is the field-free Lagrange function. Since the particle moves in 
complex time, t G [is, to], its velocity is imaginary, whereas the coordinate is real. At 
the exit x(to) = b all quantities become real. To solve f|T8l) . one has to consider t^ as an 
external parameter and tq can be found from the initial condition xo{ts = to + ^^o) = 
(for an example see the next section). 
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After the exit, when time becomes real t > to, the particle moves under the action 
of the laser field. The respective trajectory satisfies another Newton equation 

x = v = pF(t), v(to)=po, x(to) = b (20) 

with the solution 

i 

x{t) = b + p{t- to) + G{t) - G{to) , G{t) = J pF{t')dt' . (21) 



The condition x(to) = Po then specifies to fo be a function of the final momentum, p. 
The following algebra 



+ 00 , 

So = - J U/2-E0 + 



to 
+00 



d d \ . 

—vx -vx dt 

dt dt / 



(i''^/2 + vx + Eoj dt - vx\t^+oo + vx\t=to , 



to 



allows to represent the action Sq in a form identical to that of (fT9|) . Thus, the exponential 
in (|T2l) can be rewritten as exp{iW) where 
+00 

W = [ {C + Eo)dt - vx\t^+^ (22) 



ts 



is the reduced action calculated along the classical complex trajectory selected described 
above. This is the basic result of the ITM [T8] . 

The ITM provides a way to generalize the amplitude (ITTl) . Indeed, one can 
calculate the function fl22l) accounting for both the binding potential and the field of 
the electromagnetic wave, i.e. evaluating the trajectory x(t) from the equation 

x = v = -dU/dx + pF{t) (23) 

with initial and boundary conditions 

x(t,) = 0, v\ts)/2 = Eo, v{t^oo)=p. (24) 

Except for the simplest model systems, a solution to Eqs. (1231124^ can only be found 
numerically or by iteration with respect to one of the two fields. However, even in first- 
order perturbation theory, it is possible to account for the nonlinear effect of an intense 
laser field on tunneling. Indeed, if in some part of the position space the laser field is 
small compared with the binding force (or vice versa), corrections to trajectory Sx and 
to the action 5W can be derived perturbatively. These corrections must remain small 
compared to the respective unperturbed values, \6W\ -C \W\, to justify the application 
of perturbation theory. However, since under semiclassical conditions ([3]) the action is 
numerically large, the condition 

1 < \6W\ < \W\ (25) 

is typically fulfilled. This means that even a perturbative correction due to the presence 
of another (e.g. laser) field can cause a substantial modification of the spectra. The 
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regime where such a semiclassical perturbation theory for the action is relevant results 
in a variety of strong field problems [T2| [TH [15] . 

The correction SW consists of two parts: one due to the functionally different action 
that accounts for the additional interactions, and the other related to the modification 
of the trajectory. In the literature the first correction has been better studied than the 
second. In particular, in the work of Ivlev and Melnikov [19] presenting the first ITM 
treatment of laser-assisted tunneling, only the first correction was accounted for. 

To summarize, by taking into account both potential and laser fields, the differential 
probability of observing the electron with the asymptotic momentum p takes the form 
( IT71) with the amplitude 



with W{p,tsa) calculated along Eqs. f l22|) - flM|) . This distribution is the main result of 
the present work. It includes both the field-free and the laser-assisted tunneling and 
accounts for the redistribution of the particle momenta due to the laser field after exiting 
the barrier. It is relevant under conditions ([3]) with the additional requirement that the 
number of ATI-like peaks in the spectrum should be large to keep Eqs. f lTSj) and (fT6l) 
valid, i.e. PoPf oj ot zp ^ 1- Integration over the final momenta of the particle at 
the detector provides the total tunneling probability. In Appendix A it is shown that 
the field free decay rate follows from ( !26l) in the weak field limit Sq — > 0. 

A remaining important question is what happens to the spectrum for increasingly 
thick barriers up to the limit of a true bound state. For an infinitely thick barrier, the 
field-free decay vanishes and only ATI is possible. An examination of the equations of 
motion and the action considered above shows that those trajectories that correspond to 
the final energy p'^/2 > Uq survive for the infinite barrier and the respective ionization 
probability is nonzero. This will result in a common ATI spectrum (see also next 
Section). We can therefore state that the present formulation contains contributions 
from both the laser-assisted and laser-induced processes. These contributions can be 
distinguished according to the type of trajectories: trajectories that vanish for an 
infinitely thick barrier are responsible for LAT. One should note, however, that this 
classification is only qualitative, since both families of trajectories depend continuously 
on the barrier width. 

3. Rectangular barrier 

Here we consider an example of laser-assisted tunneling which admits an analytical 
solution. Let us assume a QS state in a ID rectangular well formed by short range 
forces, as depicted in Figure [TJ The central part of the well (where the particle remains 
trapped most of the time) is assumed to be narrow enough so that the external laser 
field cannot perform considerable work on the particle on this spatial scale, SqU <^ Eq. 



M{p) = 



a 



Voiyitoa)) exp(iiy(p, tsa)) 

J dp/ dtoa + i/3 



(26) 
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For simplicity we consider the limit K^a <^ 1 such that the width of the weHjjJ does not 
enter the result. 

According to the ITM recipe, one has to find trajectories that start at t = inside 
the well x{ts) = 0, with the initial kinetic energy v'^{ts)/2 = Eq — Uq = — Kq/2 (so 
that v{ts) = ino) and arrive at the detector with the particle having the desirable 
final momentum p. The detailed derivation is presented in Appendix B. For the 
monochromatic field f|T^ all periods are equivalent, so that there are only two essentially 
different solutions, all the others being obtained by a 27rj translation. The summation 
over all periods gives the factor uj5{p^/2 — Eq + Up — ju), which expresses energy 
conservation. The initial imaginary time uts = (ps = (j^o + iipo satisfies the equations 
(see also Appendix B): 

b = {kq/uj)iPq — (pp/u) cos </)o (cosh ^/'o — 1 — sinh^/^o) ? 

P = MVs)+ Pf sin 00 • (27) 

The solution of these equations and the respective trajectories given by (B.3) and 
(B.4) enter the expression of the action W{p,tsa)- The differential decay rate R over 
the whole observation time T is given by 

2 

dw{p) 2S{p-Pj) Sk-^po 



^ exp{iW{p,tsa)) 



a=l,2 



^dp/dtoa + i(3 



dp, (28) 



where Pj = yp^— p'jp/2 + 2ju are the momenta corresponding to the ATI-like peaks. 
Here we have taken into account that the prefactor Vq corresponding to the field-free 
rate Rq = PQexp(— 2Kofo) is given (in the narrow well limit kqq <^ 1) by 

4+ Pi ^ ^ 

In the weak field regime the result for the action can be simplified up to a short 
analytic form by keeping only linear terms in Sq. For a trajectory a, this procedure 
gives 

yx. / X ■ i^f, , Sobcos(j)oa , J)o^o^oosin0oa\ , Po{p-po), 

Wa{p) = tKobi 1 + + + 0o„ 

V 2kq iKQbu^ J (jj 

^0^)^ sin 00a h ■ A ^ '^Pp M ^Qn^ 

+ V'OO ^ Pb + Ppb Sm 00a H cos 0Oa , (30) 

2kq OJ 
with -00 ~ "000 = boj/nQ <^ 1. From this expression one derives the true criterion for the 
weak-field limit corresponding to the "exclusive" regime, 

/i = Kob^ < 1 . (31) 

Kq 

Under the further simplification that only the most probable trajectory is considered 
(i.e. 00 = or vr) and the real part of the action is disregarded (no interference), the 
result of Ivlev and Melnikov is recovered [19]. In this case, however, the result of [19] is 



II This should not be confused to the width of the barrier & — a « &, which is semiclassicaUy large, 
Kob^ 1. 
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more general than the one derived in this section, since the former does not assume the 
potential to be a rectangular barrier. 

For the action (l30i) . the probability vanishes in the limit 6 — )■ oo. This, however, 
cannot be consistent with the existence of ATI for an infinitely thick barrier. This 
apparent paradox can be solved if we take into account that the results of this subsection 
(including also Appendix B) are derived assuming that the particle escapes from under 
the barrier at its right edge, x = b. This is, however, true only when Sob < Uq — Eq (see 
Appendix B). Clearly, with increasing barrier thickness b, this condition will be violated 
for any given field amplitude and trajectories that escape through the tilted part of the 
barrier will come into play. Along such trajectories, the action becomes 6-independent 
and virtually identical to the case of common ATI. 

3.1. Numerical results 

As numerical examples we first consider tunneling through two rectangular barriers of 
parameters Uq = 3.0, b = 3.0 and Uq = 4.0, b = 10.0 assisted by a laser field with 
frequency u = 0.1. The action W and the rate R are obtained using the expressions 
given in Appendix B and Eq. fl28|l . In Figure [2] we present the imaginary part of 
the action and the spectrum of the laser-assisted tunneling as a function of the final 
energy E for three different field amplitudes Sq. For each energy between -Emin and 
-^'max in Figure [2](a) and (d) the imaginary part of the action Im(W) has two values 
corresponding to the two trajectories starting inside the well at times ipsa with a = 1,2. 
The spectrum consists of several ATI-like maxima whose positions are dictated by the 
energy conservation conditions in Eq. ( l28l) . 

In the case of the thin barrier with 6 = 3, the tunneling occurs mostly field-free, 
and the main effect of the laser is to change the particle momentum after the barrier 
exit. By comparing the results for different field intensities in Figure [2)d, we observe that 
the spectra become more narrow with decrease of Sq, approaching the Lorenzian shape 
of the field-free decay rate. Furthermore, the ratios between the field-free and laser- 
induced tunneling probabilities do not vary much for the considered field amplitudes 
Rom/Ro = 0.41, Rq^q^/Rq = 0.62 and Rq.u/Rq = 0.89. Note that for these parameters, 
the large number of ATI-like maxima condition poPr ^ i^, necessary to justify the 
momentum distribution along Eq. (fTSjl is not fulfilled. As a consequence, our distribution 
loses a part of the tunneled particles, and the laser-assisted to field-free decay ratios are 
less than unity. Due to this effect, it is more informative to consider as reference the 
laser-assisted rate at the lowest field intensity. This gives us -R0.12/-R0.02 ~ 2, showing 
that the laser field of the amplitude Sq = 0.12 only enhances the total decay rate by the 
factor of 2 for this barrier width. 

In contrast, for the thick barrier of width b = 10, corresponding to a very small field- 
free decay rate, the "inclusive" regime is achieved, in which the laser has a substantial 
effect also on the tunneling rate itself. Here, the ratios of the field-free and laser-induced 
tunneling probabilities are -R0.02/-R0 = 0.72, Rqq^/Rq = 1.94 and Rq ^/Ro = 58.17, 
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Figure 2. Imaginary part of the action Ini(W^) (a), (c) and the spectra (b), (d) as a 
function of the the final energy E = for two sets of parameters: Uq = 3.0, b ~ 3.0, 
Eo = 1.217, CO = 0.10, (a)-(b) and Uq = 4.0, b = 10.0, Eq = 1.302, tj = 0.10, (c)-(d). 
Each panel shows three curves for Sq = 0.02 (black double-dashed line). So = 0.05 
(dashed green line) and £q — 0.12 (solid red line). 

extending over almost two orders of magnitude with the increasing field. Thus, 
as anticipated in the introduction, our approach reproduces important qualitative 
conclusions of earlier studies, namely that depending on the parameters, there are two 
different regimes of decay: "exclusive" , when the spectrum is strongly affected without 
a modification of the total decay rate (the thin barrier) and "inclusive" , when the rate of 
decay is strongly modified by the laser field (the thick barrier). Quantitatively, these two 
regimes can be distinguished according to value of the parameter ( ISTl) . The "exclusive" 
regime requires yU -C 1 and holds for the first set of parameters, while for the thick 
barrier we enter the "inclusive" domain, ^ 1. 

In Figure |3] we have considered the case of the thin barrier with b = 3 for different 
laser field amplitudes and frequencies such that the field momenta are pp = 1. The ratio 
of the different field frequencies can be identified from the spectra. The total tunneling 
rate is increasing with increasing the field, with Rq i/Rq = 0.74, Rq2/Ro = 0.91 and 
Ro.s/Ro = 1.56. 

If we consider instead of the monochromatic field a few-cycle laser pulse, the 
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Figure 3. Imaginary part of the action (left) and laser assisted decay spectra (right) 
calculated for monochroniatic fields with various frequencies and field momenta pj;^ = 1 
for the thin barrier considered in Figure [2l The field amplitudes and frequencies are 
So = 0.1 and w = 0.1 (black double-dashed line), Sq = 0-2 and uj = 0.2 (green dashed 
line) and So = 0.3 and uj = 0.3 (red solid line), respectively. 



spectrum of Dirac-delta functions will be replaced by a comb of broadened ATI-like 
maxima between E'min and i^^max- Spectra obtained using the extended ITM for a six- 
cycle pulse of the form 

= £^0 sin^ I J coswt , Up = Q (32) 

for the parameters previously addressed in the text are presented in Figure HI Here the 
amplitude consists of up to 2np coherent contributions that produce the interference 
pattern of the spectra. The broad ATI-like maxima can be observed best for the thick 
barrier case in Figure HJi, where the absolute values of the two contributions from 
a given laser period differ substantially and smear out the interference. Unlike the 
case of tunneling assisted by a monochromatic field where the differential rate dR{p) is 
calculated, the spectra in Figure H] present the differential probability dw{p)/dp. The 
ratio w(p) /pulse duration delivers approximately the decay rate. 

3.2. Comparison with numerical results of the time- dependent Schrddinger equation 
(TDSE) 

In order to check the accuracy of our approach we compared the obtained spectra with 
accurate numerical results of the TDSE. For the test case of the rectangular barrier, we 
have calculated exact numerical spectra using a ID version of the QPROP Schrodinger 
solver [33], which propagates the wavefunction in real time on a spatial grid. For the 
numerical simulations, the width of the potential well (in which f/ = 0) needed to be 
specified and was chosen as a = 7r/2. We first obtained the ground state wavefunction 
for the particle inside the well of height f/o=3 and h — )■ oo, i.e. a barrier of infinite width, 
as depicted by the dashed line in Figure [H The ground state energy on the numerical 
grid of spacing Ax = 0.1 was Eq = 1.24. At time t = 0, the infinite barrier was replaced 
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Figure 4. Imaginary part of tire action (a), (c) and laser assisted decay spectra (b), 
(d) calculated for a finite pulse of the form (|32p for the two barrier parameter sets of 
Figure [2] and fields of frequency uj — 0.1 and amplitudes £o = 0.02, 0.05, 0.12. 

by a barrier of finite width b — a = 4, keeping tlie height constant. As a consequence, 
tunnehng occurs for t > 0. The sudden switch from the infinite to the finite barrier 
disturbs the trapped electron and leads to a short time interval of transient tunneling 
dynamics before a constant free-tunneling rate is established. However, this time interval 
of a few atomic units is much shorter than the pulse duration so that it did not make 
a difference whether the six-cycle laser pulse of the form f l32|) and frequency u = 0.057 
(Ti:Sa laser) was switched on at t = or with a delay. In any case, the numerical grid 
was big enough to support the entire wavefunction during the propagation time without 
reflections off the grid boundary. The electron spectra were calculated using the window 
operator technique (see, e.g., [33]) and normalized to the field- free decay. 

In Figure O we compare the obtained tunneling probabilities for a barrier of 
thickness b — a = 4, height Uq = 3 and initial particle energy Eq = 1.24 under the 
action of fields of amplitudes £o=0.02 (Figure |5t) and £^o=0.005 (Figure |5b). The 
spectra agree well both qualitatively and quantitatively for final energies within the 
CB, with the ITM results slightly higher than the TDSE ones. From a comparison of 
the two field-free decay rates, we observe that the fitted from the TDSE results are 
always smaller than the calculated Rq. This behavior is related to TDSE numerics 
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Figure 5. (a),(b),(c) Tunneling probabilities dw{E)/dE for rectangular barriers of 
widths b — a = 4 (a), (b), and 5 — a = 6 (c) as a function of the final energy: the 
present ITM calculation (solid red line) and numerical solution of the TDSE (green 
dashed line). The laser field parameters are uj — 0.057 and (a) £o=0.02, (b) £f)~0.005 
and (c) £o~0.05. (d) Comparison between laser-assisted tunneling through the finite 
barrier considered in (a) under the action of a field of amplitude £o=0.075 (solid red 
line) and ATI for 6 = oo in the same field (green dashed line). 



requirements, which cannot handle the very thin potential well limit a — )■ 0. 

The ITM approach delivers spectra that vanish abruptly beyond the CB and cannot 
reproduce the shoulders visible in the TDSE results. As already discussed in Section |2TT| 
the saddle-point method is actually not applicable in the form described here outside 
the CB. The correct approach requires us to include the term proportional to the third 
derivative of the action in the phase decomposition near the saddle point [32]. For 
broad spectra, the contribution of the shoulders outside the CB is not significant and 
our approach provides reliable results. The contribution of the shoulders increases for 
narrow spectra, as one can see comparing Figures [SH and [5b- 

We further compare the ITM and TDSE results for a thicker barrier with 6 — a = 6 
and field amplitude £^o=0.05 in Figure |5t. Here the agreement is less accurate at small 
energies E < 0.8, where the ITM results are about one order of magnitude higher than 
the TDSE ones. Since ITM delivers a momentum spectrum, the variable transformation 
dE = pdp introduces a divergence for asymptotic momenta approaching the origin. 
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However, for energies close to the initial particle energy and larger, the ITM agrees 
well with the TDSE. We conclude that the ITM provides not only qualitative but also 
quantitative results for laser-assisted tunneling of QS states within the semiclassical 
parameter regime. 

Note that in all numerical examples considered in this section, two effects were 
disregarded in the ITM. Firstly, in our analytical formulae trajectories that escape 
through the tilted barrier are not accounted for, so that the contribution that turns 
in the limit 6 — )■ oo into ATI is missing. For the parameters we have chosen, this 
contribution can be safely neglected. For demonstration, we show in Fig. [Sji the LAT 
spectrum calculated numerically using the TDSE for the barrier parameters of panel 
(a) under the action of a laser field of amplitude Sq = 0.075 and compare it with ATI 
through an infinite barrier {b = oo) of the same height. We can see that for the chosen 
parameters, the ATI probability is many orders of magnitude smaller than laser-assisted 
tunneling one. Secondly, the theory does not take into account that the electron can be 
driven back to the barrier and rescatter absorbing or emitting additional photons. As is 
known from the literature, rescattering leads to the formation of one or more plateaus 
in the spectrum, with the characteristic number of peaks given by the Reiss parameter 
(jl]). The same effect can also be interpreted as multiphoton stimulated bremsstrahlung 
[2]. For our calculations, we have selected the parameters such that in all cases Pq > pp 
and rescattering plays no significant role. The rescattering plateau can be reproduced 
by TDSE calculations covering the higher field amplitude domain pp > pq, which was, 
however, not considered in this work. 

4. Conclusions 

Starting from the SFA and its formulation in terms of complex trajectories, we have 
developed a general method to describe the laser-assisted decay of QS states in 
the semiclassical parameter regime. In order to illuminate the physical essence of 
the problem and avoid unnecessary technical complications, we have neglected some 
accompanying effects like conventional ATI and rescattering. A test case comparison 
with numerical results of the Schrodinger equation shows not only qualitative but also 
quant it at i ve agreement . 

Due to the general statement of the problem, our method can be applied 
with moderate numerical effort to many laser-assisted tunneling processes within the 
semiclassical regime. In particular, laser-assisted tunneling through a more realistic 
Coulomb barrier as it occurs in autoionization, autoemission or generation of high 
harmonics in the presence of static or low-frequency fields [34] can also be described 
via our extended ITM. Here the real strength of the ITM could be used to correctly 
take into account the effect of long-range potentials that are acting on the particle in the 
classically allowed domain simultaneously with the laser field. Work in this direction is 
in progress. 

As a further incentive, recently ATI was experimentally observed at sharp metal 
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tips driven by moderately intense infrared laser pulses in the presence of a DC field [55] . 
The latter formed an effective tilted barrier that the electrons tunneled through with 
the assistance of the laser field. Under the action of the AC laser field, the electron 
current was modulated in time, forming femto- and attosecond dense electron bunches 
[36] . An extension beyond the ID case of the method developed in this work can be 
applied to model such experiments. For non-separable variables in many dimensions, 
our approach is expected, however, to involve more cumbersome calculations (see also 
Ref. [37]), since the corresponding multidimensional trajectories and complex times can 
only be found numerically. This should, however, remain an essentially simpler task as 
compared to obtaining the numerical solution of the corresponding TDSE. 



Acknowledgements 

We would like to thank C. H. Keitel, W. Becker, V. D. Mur and V. S. Popov for fruitful 
discussions. AP acknowledges financial support of the Robert-Bosch-Stiftung in the 
framework of the "Fast Track" Programme. SVP acknowledges financial support from 
the Deutsche Forschungsgemeinschaft (project No. SM 292/1-1) and from the Russian 
Foundation for Basic Research (project No. 11-02-91331). DB acknowledges support 
from the Deutsche Forschungsgemeinschaft (SFB 652). 



Appendix A. Field-free limit 

In this Appendix we show that in the weak field limit the field-free decay rate 

i?o = Vq exp(-2Imiyo) (A.l) 
follows from the amplitude ( l26l) . Quantitatively, the weak field limit is determined by 
the condition fl3T]) . Correspondingly, the subbarrier correction to the action is smaller 
than unity so that the imaginary part of the action is given by the field-free contribution 
Wq ( !T9|) . Since the momentum change during the subbarrier motion is also small, the 
initial velocity at the exit is v{to) = po and then the final momentum is given by 

Pito) = po + Pf sin ut . (A.2) 

We have assumed here a monochromatic field (HM for simplicity. The spectrum consists 
of L = 2pqPp/u ^ 1 ATI-like peaks with energies between (po+Pf)^/2 and {po—pp)^/2. 
For the given momentum inside this interval, there are two solutions per period, so that 

^ton = arcsin ( - — — ] + 27rn , 
\ Pf J 

ooto^^ = TT — a.Tcsm(- — —)+27rn. (A. 3) 



Pf 

In the limit we consider, the laser field only enters the action via these initial times tQ, 

T 



W = Wo + J (pV2 + Eo)dt - px{T) + poh 



to 



ilmWo + ^{pl-p'){T-to), (A.4) 
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where T is the large observation time and we take into account that x{t) ^ 6 + — tg) 
The sum over the laser periods gives 



N 



Then 



n=0 



dR 



^exp [i{p'^ -pI 
VIM 



Tin 



3 



N ^ oo . 



(A.5) 



dw 



exp(-2Imiyo 



X 



vl 



2uj 



J 



exp{ijujtl^^] 



exp(-zjw4 



(+)^ 



COS cut 







— COS ut. 







dp. (A.6) 



Under the condition p^pp/u ^ 1, the number of ATI-like peaks is large and the sum 
over j can be replaced by an integral which evaluates to 1. The resulting distribution 
should be integrated over dp within the limits po ± Pf- Taking into account that 
cos utQ = — {p — PqY Ip^p and disregarding the rapidly oscillating interference term 
in the modulus square, we obtain precisely the field- free rate (lA.ll) . 



Appendix B. Rectangular barrier 

Here we present in more detail the ITM treatment for the case of laser-assisted tunneling 
through a rectangular barrier discussed in Section [31 In the field-free case, to each initial 
time instant to corresponds the trajectory 

Xo{(p) = b+ {po/u){(p - (po) , V5e[0o,+oo)- (B.l) 

Here and below we use the dimensionless time if = ut = (p + iip such that ips = 0o + ^"^0- 
The field-free action is then 

Wo = ibKo-pob. (B.2) 

The trajectory is always real, while the velocity is imaginary in complex time and 
becomes real when time is real at x > 6. This allows for the standard interpretation that 
the particle moves under the potential barrier having an imaginary velocity. Later on, 
the particle escapes from under the barrier and its trajectory becomes entirely classical. 
As it will be seen in the following, for the time-dependent case this subdivision no longer 
holds: a trajectory can have a nonzero imaginary part even in real time. However, for 
the sake of convenience, we will still refer to subbarrier motion when time is complex, 
ip ^ and to motion after the barrier otherwise. 

In the presence of the field (JT^ . the trajectory under the barrier is given by 

v\(p) = Ko[i - 7"^(sinv9 - sinip^)] , 

x\(p) = -^[{h + smips){(p - ips) + (cosv5 - cosips)] . (B.3) 
After the barrier, if > (po, where U{x) = 

f ™(</)) = vq — ppisin ip — sin 0o) , with p = Vq + pp sin 0o , 

= x^(0o) +p/^{^ - 0o) + (pf/w)(cosv? - COS0o) • (B.4) 
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Here we introduce the dimensionless parameter 

7 = kquj/So (B.5) 

having the same physical meaning as the Keldysh parameter [5], i.e. the ratio of the 
"atomic" momentum to the field- induced momentum pp = Sq/u. 

Trajectories (IB.Sp and (]B.4p are, however, insufficient to find the initial complex 
time, because the velocity is discontinuous at the barrier exit, x = b. To avoid this 
discontinuity we replace the step-like potential by a smooth one; for example, the one 
depicted in Figure [T] by the dashed-dotted line. The potential then drops from Uq 
linearly down to zero on the width A. This triangular barrier we denote as region II, 
while outside the barrier we define region III. In the intermediate region II the solution 
is 

v^\ip) =v\if') -pFisimp-smip') + {Fo/u){ip - ip') , Fq = f/o/A , 

= x\<p') + {pf/uj){cos<p - cos<p') + {Fo/2u^){<p - <p'f . (B.6) 

Here (f' is the time instant when Re[x^ = b. Matching the solutions in the domains 
I, II and HI, we can find the constant vq as a function of the other parameters. Note that 
the sharper the slope is, the closer the time instant ip' = (f)o + iip' is to 00 . Decomposing 
the equations with respect to ip' we obtain: 

■^0 = (Po + 2koPf cos 00 sinh 

1 /2 

+ p|'[(cosh'?/'o — 1)^ sin^ 00 — cos^ 00 sinh^ ■i/'o] f • (B-7) 



The two equations determining the initial complex time follow from the requirements 
x(0o) = b and w(0 — ?■ oo) = p and take the form ( 127|) . They can be solved analytically in 
the weak field limit, Sq/oj'^ <^ or low-frequency field 7 <^ 1, while in the general case 
the numerical solution is required. Furthermore, solutions do not exist for all values of 
the final momentum and arbitrary parameters and 7. Indeed, the second equation 
of ( !27l) shows that the final momentum is determined by the time instant 0o when the 
particle is released from under the barrier, and by its initial velocity vq. The width of 
the momentum space available is determined by the field momentum, Pmax.min ~ Vo^Vfi 
so that when po — Pf < not all 0o are solutions. Another limitation comes from the 
fact that Eqs. (B.3)-(B.6) assume that the instant kinetic energy at the exit is below 
the instant barrier height, Vq < 2{Uq — Sob cos (po); this restricts the value of the field. 
Sob < K,l/2. 

In the rectangular barrier limit A — 0, the domain II vanishes and the action is 
determined by trajectories ( IB. 31) and (IB. 41) and the Lagrange function, 

C = xV2 + SoX cos (fi — UoQib — x) , (B.8) 

where G(x) is the Heaviside function. Using this Lagrangian and the above determined 
trajectories, the classical action has to be evaluated along Eq. (122|) . 



Laser assisted decay of quasistationary states 



24 



References 

[1] Voronov G S, Delone N B 1965 JETP Lett. 1 66 (Engl. Transl.) 

[2] Bunkin V F, Fedorov M V 1966 Sov. Phys. - JETP 22 844 (Engl. Transl.) 

[3] Fedorov M V 1997 Atomic and Free Electrons in a Strong Ligtit Field (Singapore: World Scientific) 
[4] Manakov N L, Frolov M V, Starace A F and Fabrikant 1 1 2000 J. Phys. B: At. Mol. Opt. Ptiys. 33 
R141 

[5] Keldysh L V 1965 Sov. Pfiys. JETP 20 1307 [1964 Zh. Eksp. Teor. Fiz. 47 1945] 
[6] Faisal F H M 1973 J. Ptiys. B 6 L89 
[7] Reiss H R 1980 Phys. Rev. A 22 1786 
[8] Popov V S 2004 Phys.-Usp. 47 855 

[9] Milosevic D B, Paulus G G, Bauer D and Becker W 2006 J. Phys. B: At. Mol. Opt. Phys. 39 R203 
[10] Volkov D M 1935 Z. Phys. 94 250 

[11] Berestetsky V B, Lifshits E M, Pitaevsky L P 1971 Quantum electrodynamics (Oxford: 

Butterworth-Heinemann) 
[12] Perelomov A M and Popov V S 1967 Sov. Phys. JETP 25 336 (Engl. Transl.) 
[13] Popov V S, Kuznetsov V P, Perelomov A M 1967 Sov. Phys. JETP 53 331 (Engl. Transl.) 
[14] Popruzhenko S V, Mur V D, Popov V S and Bauer D 2008 Phys. Rev. Lett. 101 193003; JETP 

2009 108 947 [2009 Zh. Eksp. Teor. Fiz. 135 1092] 
[15] Popruzhenko S V, Paulus G G, Bauer D 2008 Phys. Rev. A 77 053409; Popruzhenko S V, Bauer 

D 2008 J. Mod. Opt. 55 2573 
[16] Huismans Y et al. 2011 Science 331 61 

[17] Yan T-M, Popruzhenko S V, Vrakking M J J and Bauer D 2010 Phys. Rev. Lett. 105 253002 
[18] Popov V S 2005 Phys. At. Nuclei 68 686 

[19] Ivlev B I, Melnikov V I 1985 JETP Letters 41 116 (Engl. Transl.) 
[20] Platero G, Aguado R 2004 Phys. Rep. 395 1 

[21] Nikishov A I, Ritus V I 1964 Sov. Phys. JETP 19 1191 [1964 Zh. Eksp. Teor. Fiz. 46 1768] 

[22] Becker W, Moore G T, Schhcher R R and Scully M O 1983 Phys Lett A 94 131 

[23] Becker W, Schhcher R R and Scully M O 1984 Phys Lett A 101 58 

[24] Becker W, Schhcher R R and Scully M O 1984 New Phys A 426 125 

[25] Fedotov A M, Narozhny N B, Mourou G, and Korn G 2010 Phys. Rev. Lett. 105 080402 

[26] Bauer D, Milosevic D B, Becker W 2005 Phys. Rev. A 72 023415 

[27] Smirnova O, Spanner M, Ivanov M 2008 Phys. Rev. A 77 033407 

[28] Baz A I, Zeldovich Ya B, Perelomov A M 1969 Scattering, Reactions and Decays in Nonrelativistic 

Quantum Mechanics (Jerusalem: Israel Program for Scientific Translations Press) 109-46. 
[29] Peierls R 1991 More surprises of quantum mechanics (Princeton: University Press) 
[30] Gamow G 1928 Z Phys 51 204 

[31] Zeldovich Ya B 1961 Sov. Phys.-JETP 12 542 (Engl. Transl.) 

[32] Goreslavski S P and Popruzhenko S V 1999 J. Phys. B: At. Mol. Opt. Phys. 32 L531; Figueira de 

Morisson Faria C, Schomerus H and Becker W 2002 Phys. Rev. A 66 043413 
[33] Bauer D, Koval P 2006 Comput. Phys. Commun. 174 396 

[34] Bao M-Q and Starace A F 1996 Phys. Rev. A 53 R3723; Odzak S and Milosevic D B 2005 Phys. 
Rev. A 72 033407 

[35] Schenk M, Kriiger M and Hommelhoff P 2010 Phys. Rev. Lett. 105 257601 

[36] Papadichev V A 2011 private communication 

[37] Kapur P R, Peierls R 1937, Proc. Roy. Soc. A163 606 



